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Severe acute respiratory syndrome (SARS) 
was caused by a newly emerged coronavims, 
now known as SARS coronavims (SARS-CoV) 
(7, 2). In spite of the early success of etio¬ 
logical studies and molecular characterization 
of this vims (3, 4 ), efforts to identify the ori¬ 
gin of SARS-CoV have been less successful. 
Without knowledge of the reservoir host dis¬ 
tribution and transmission routes of SARS- 
CoV, it will be difficult to prevent and control 
future outbreaks of SARS. 

Studies conducted previously on animals 
sampled from live animal markets in Guang¬ 
dong, China, indicated that masked pahn civets 
(Paguma larvata) and two other species had 
been infected by SARS-CoV (5). This led to 
a large-scale culling of civets to prevent fur¬ 
ther SARS outbreaks. However, subsequent 
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studies have revealed no widespread infection 
in wild or farmed civets (6, 7). Experimental 
infection of civets with two different human 
isolates of SARS-CoV resulted in overt clin¬ 
ical symptoms, rendering them imlikely to be 
the natural reservoir hosts (#). These data sug¬ 
gest that although P. larvata may have been 
the source of the human infection that pre¬ 
cipitated the SARS outbreak, infection in this 
and other common species in animal markets 
was more likely a reflection of an “artificial” 
market cycle in naive species than an indi¬ 
cation of the natural reservoir of the virus. 

Bats are reservoir hosts of several zoonotic 
viruses, including the Hendra and Nipali viruses, 
which have recently emerged in Australia and 
East Asia, respectively (9-11). Bats may be 
persistently infected with many viruses but rarely 
display clinical symptoms (12). These charac¬ 
teristics and the increasing presence of bats and 
bat products in food and traditional medicine 
markets in southern China and elsewhere in 
Asia (13) led us to survey bats in the search 
for the natural reservoir of SARS-CoV. 

In this study, conducted from March to 
December of 2004, we sampled 408 bats 
representing nine species, six genera, and 
three families from four locations in China 
(Guangdong, Guangxi, Hubei, and Tianjin) af¬ 
ter trapping them in their native habitat (Table 
1). Blood, fecal, and throat swabs were col¬ 
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lected; serum samples and cDNA from fecal 
or throat samples were independently analyzed, 
double-blind, with different methods in our 
laboratories in Wuhan and Geelong (14). 

Among six genera of bat species surveyed 
(Rousettus, Cynopterus, Myotis, Rhinolophus, 
Nyctalus, and Miniopterus), three communal, 
cave-dwelling species from the genus Rhinolo- 
phus (horseshoe bats) in the family Rhino- 
lophidae demonstrated a high SARS-CoV 
antibody prevalence: 13 out of 46 bats (28%) 
in R. pearsoni from Guangxi, 2 out of 6 bats 
(33%) in R. pussilus from Guangxi, and 5 out 
of 7 bats (71%) in R. macrotis from Hubei. 
The high seroprevalence and wide distribution 
of seropositive bats is expected for a wildlife 
reservoir host for a pathogen (15). 

The serological findings were corroborated 
by poylmerase chain reaction (PCR) analyses 
with primer pairs derived from the nucleocapsid 
(N) and polymerase (P) genes (table SI). Five 
fecal samples tested positive, all of them from 
the genus Rhinolophus : three in R. pearsoni 
from Guangxi and one each in R. macrotis and 
R. ferrumequinum, respectively, from Hubei. 
No virus was isolated from an inoculation of 
Vero E6 cells with fecal swabs of PCR-positive 
samples. 

A complete genome sequence was deter¬ 
mined directly from PCR products from one of 
the fecal samples (sample Rp3) that contained 
relatively high levels of genetic material. The 
genome organization of this vims (Fig. 1), ten¬ 
tatively named SARS-like coronavims isolate 
Rp3 (SL-CoV Rp3), was essentially identical to 
that of SARS-CoV, with the exception of three 
regions (Fig. 1, shaded boxes). The overall 
nucleotide sequence identity between SL-CoV 
Rp3 and SARS-CoV Tor2 was 92% and in¬ 
creased to ~94% when the three variable re¬ 
gions were excluded. The variable regions are 
located at the 5" end of the S gene (equivalent to 
the SI coding region of coronavims S protein) 
and the region immediately upstream of the 
N gene. These regions have been identified 
as “high mutation” regions among different 
SARS-CoVs (5, 16, 17). The region upstream 
of the N gene is known to be prone to de¬ 
letions of various sizes (J, 16, 18). 

Predicted protein products from each gene 
or putative open reading frame (ORF) of SL- 
CoV Rp3 and SARS-CoV Tor2 were com- 
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Table 1. Detection of antibodies to SARS-CoV and PCR amplification of N and P gene fragments with 
SARS-CoV-specific primers. ND, not determined because of poor sample quality or unavailability of 
specimens from individual animals. 



Sampling 

Bat species 

Antibody test: 
positive/total (%) ■ 

PCR analysis: 
positive/total (%) 

Time 

Location 

Fecal Respiratory 

swabs swabs 

Mar 04 

Nanning, Guangxi 

Rousettus leschenaulti 

1/84 (1.2%) 

0/110 

ND 

Maoming, Guangdong 

Rousettus leschenaulti 

0/42 

0/45 

ND 



Cynopterus sphinx 

0/17 

0/27 

ND 

July 04 

Nanning, Guangxi 

Rousettus leschenaulti 

ND 

0/55 

0/55 


Tianjin 

Myotis ricketti 

ND 

0/21 

0/21 

Nov 04 

Yichang, Hubei 

Rhinolophus pusillus 

ND 

0/15 

ND 



Rhinolophus ferrumequinum 0/4 

1/8 (12.5%)* 

ND 



Rhinolophus macrotis 

5/7 (71%) 

1/8 (12.5%)f 

0/3 



Nyctalus plancyi 

0/1 

0/1 

ND 



Miniopterus schreibersi 

0/1 

0/1 

ND 



Myotis altarium 

0/1 

0/1 

ND 

Dec 04 

Nanning, Guangxi 

Rousettus leschenaulti 

1/58 (1.8) 

ND 

ND 



Rhinolophus pearsoni 

13/46 (28.3%) 

3/30 (10%)$ 

0/11 



Rhinolophus pussilus 

2/6 (33.3%) 

0/6 

0/2 

*Positive 

fecal sample designated 

Rfl fPositive fecal sam 

iple designated Rml 

JPositive fecal samples 

designated Rpl, Rp2, and Rp3, respectively. 
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Fig. 1. Genome organization of, and comparison between, SL-CoV and SARS-CoV. (A) Overall 
genome organization of SL-CoV RpB. (B) Expanded diagram of the 3' region of the genome in 
comparison with SARS-CoV strains Tor2 and SZ3, following the same nomenclature used by Marra 
et al. (4). The genes (named by letters P, S, E, M, and N) present in all coronaviruses are shown in 
dark-colored arrows, whereas the SARS-CoV-specific ORFs are numbered and illustrated in light- 
colored arrows. ORFIO' follows the nomenclature by Guan et al. (5) to indicate that the single ORF 
present between ORF9 and N in SL-CoV is equivalent to the fusion of ORFIO and ORF11 in the 
same region in SARS-CoV Tor2. The shaded boxes mark the only three regions displaying 
significant sequence difference between the two viruses (table S2). 


pared (table S2). The P, S, E, M, and N pro¬ 
teins, which are present in all coronaviruses, 
were similarly sized in the two viruses, with 
sequence identities ranging from 96% to 
100%. The only exception was the SI domain 
of the S protein, where sequence identity fell to 
64%. The SI domain is involved in receptor 
binding, whereas the S2 domain is responsible 
for the fusion of virus and host cell membranes 
(19). The sequence divergence in the SI do¬ 
main corroborated our serum neutralization 


studies, which indicated that although bat sera 
have a high level of cross-reactive antibodies 
(with enzyme-linked immunosorbent assay 
titers ranging from 1:100 to 1:6400), they 
failed to neutralize SARS-CoV when tested on 
Vero E6 cells. This finding suggests that SI is 
the main target for antibody-mediated neutral¬ 
ization of this group of viruses, which is con¬ 
sistent with previous reports indicating that 
major SARS-CoV neutralization epitopes are 
located in the SI region (20, 21). 


In addition to the five genes present in all 
coronavirus genomes, coronaviruses also have 
several ORFs between the P gene and the 3' 
end of the genome that code for nonstructural 
proteins. The function of these nonstructural 
proteins is largely unknown. The location and 
sequence of ORFs are group- or virus-specific 
and hence can serve as important molecular 
markers for studying virus evolution and clas¬ 
sification (19, 22). SARS-CoV has a unique 
set of ORFs not shared by any of the known 
coronaviruses (3, 4). Most of these ORFs 
were also present in SL-CoV, confirming the 
extremely close genetic relationship between 
SARS-CoV and SL-CoV (Fig. 1 and table S2). 

Coronaviruses produce subgenomic mRNAs 
through a discontinuous transcription process 
not fully characterized (19). Conserved nucle¬ 
otide sequences functioning as transcription 
regulatory sequences (TRSs) are required for 
the production of the subgenomic mRNAs. In 
SARS-CoV, such TRSs were identified at 
each of the predicted gene start sites (3, 4). 
All of these TRSs were absolutely conserved 
between SARS-CoV Tor2 and SL-CoV Rp3 
(table S3), further demonstrating that these 
two viruses are very closely related. 

SL-CoV is completely different from a 
bat coronavirus (bat-CoV) recently identified 
by Poon et al. (7) from species of bats in the 
genus Miniopterus during a wildlife surveil¬ 
lance study in Hong Kong (Fig. 2). Because 
the complete genome sequence was not 
available for bat-CoV, only the trees covering 
the common sequences (i.e., parts of the Plb 
and S2 proteins) are shown. The phylogenetic 
analysis demonstrated that SL-CoV Rp3 and 
SARS-CoVs are clustered together but that 
bat-CoV is placed among the relatively dis¬ 
tant group 1 viruses. Hereafter, SARS-CoVs 
and SL-CoVs will be collectively called the 
SARS cluster of coronaviruses. 

In addition to the complete genome se¬ 
quence of SL-CoV Rp3, partial genome se¬ 
quences for the other four PCR-positive bat 
samples were also determined. Phylogenetic 
analysis based on the N protein sequences 
(Fig. 3A) revealed that the genetic variation 
among the SL-CoV sequences was much greater 
than that exhibited by SARS-CoVs (for sim¬ 
plicity, only three human and civet SARS-CoV 
isolates were used; the remainder are ahnost 
identical to those shown). This was especially 
obvious when SL-CoVs isolated from different 
bat species were compared. Moreover, the re¬ 
sults suggested that SARS-CoVs nestle phylo- 
genetically within the spectrum of SL-CoVs. 

We also compared the “high mutation” re¬ 
gions in samples Rfl, Rml, and Rp3. For the 
region upstream of the N gene, SL-CoVs from 
all three bat species contained a single ORF 
(ORFIO'), similar to that found in SARS-CoV 
isolates from civets (5) and patients in the 
early phase of the outbreaks (16, 18) but dif¬ 
ferent from that in most human isolates, which 
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Fig. 2. Phylogenetic 
trees. (A) and (B) are 
trees based on de¬ 
duced amino acid se¬ 
quences of the same 
regions in PIb and S, 
respectively, as used by 
Poon et al. (7) for bat- 
CoV. Tor2 and SZ3, 
SARS-CoV strains Tor2 
and SZ3; Rp3, SL-CoV 
Rp3; HCoV, human co- 
ronavirus; MHV, mouse 
hepatitis virus; PEDV, 
porcine epidemic diar¬ 
rhea virus; IBV, avian in¬ 
fectious bronchitis virus. 
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Fig. 3. Phylogenetic trees based on deduced ami¬ 
no acid sequences of (A) N, (B) ORF10', and (C) 
SI proteins. Tor2, SZ3, and GD01, different SARS- 
CoV strains; Rfl, Rml, and Rpl-3, different SL- 
CoV sequences. The genetic distance scale shown 
for (A) is different from those for (B) and (C). 


have a 29-nucleotide deletion in this region 
(3, 4, 16). ORF10' in Rfl codes for a protein 
having the same size (122 amino acids) as and 
more than 80% sequence identity to ORF10' 
proteins of SARS-CoVs, but those in Rml 
and Rp3 code for a 121-amino acid protein 
with only 35% sequence identity (Fig. 3B and 
fig. S2). By contrast, analysis of the SI pro¬ 
tein regions (Fig. 3C and fig. S3) indicated 
that Rfl was more closely related to SL-CoVs 
from two other bat species than to SARS- 
CoVs, suggesting that the SARS cluster of 
coronaviruses could recombine to increase 
genetic diversity and fitness, as is well doc¬ 
umented for other coronavimses (19). We 
were unable to sequence these regions for Rpl 
or Rp2, owing to the poor quality of the fecal 


materials from these two animals. The limited 
amount of cDNA available was used up for N 
gene analysis and in initial sequencing trials 
with SARS-CoV-derived primers, which 
were largely unsuccessful. Judging from the 
close relationship of the N genes between 
Rpl, Rp2, and Rp3 (fig. SI), it is unlikely that 
Rpl or Rp2 will have major sequence differ¬ 
ences from Rp3 in the SI or ORFIO' regions. 
This is not unexpected, considering that these 
three positive samples were obtained from the 
same bat species in the same location. 

The genetic diversity of bat-derived se¬ 
quences supports the notion that bats are a 
natural reservoir host of the SARS cluster of 
coronaviruses. A similar observation has been 
made for henipaviruses, another important group 
of emerging zoonotic viruses of bat origin, 
which show greater genetic diversity in bats 
than was observed among viruses isolated 
during the initial Nipah outbreaks in Malaysia 
(23-26). The overall nucleotide sequence iden¬ 
tity of 92% between SL-CoVs and SARS- 
CoVs is very similar to that observed between 
Nipah viruses isolated from Malaysia and 
Bangladesh in 1999 and 2004, respectively 
(25) (fig. S4). SL-CoVs present a new chal¬ 
lenge to the diagnosis and treatment of future 
disease outbreaks. The current tests and ther¬ 
apeutic strategies may not work effectively 
against all viruses in this group, owing to 
their great genetic variability in the SI do¬ 
main region of the S gene. 

The genus Rhinolophus contains 69 species 
and has a wide distribution from Australia to 
Europe (27). They roost primarily in caves and 
feed mainly on moths and beetles. However, 
notwithstanding the predominant Rhinolophus 
findings in this study, it is highly likely that 
there are more SARS-related coronaviruses to 
be discovered in bats. Indeed, our positive 
serological findings in the cave-dwelling fruit 
bat Rousettus leschenaulti indicate that infec¬ 


tion by a related virus could occur in fruit bats 
as well, albeit at a much lower frequency. A 
plausible mechanism for emergence from a 
natural bat reservoir can be readily envisaged. 
Fruit bats including R. leschenaulti, and less 
frequently insectivorous bats, are found in mar¬ 
kets in southern China. An infectious consign¬ 
ment of bats serendipitously juxtaposed with a 
susceptible amplifying species, such as P. 
larvata, at some point in the wildlife supply 
chain could result in spillover and establish¬ 
ment of a market cycle while susceptible ani¬ 
mals are available to maintain infection. Further 
studies in field epidemiology, laboratory infec¬ 
tion, and receptor distribution and usage are 
being conducted to assess potential roles played 
by different bat species in SARS emergence. 

These findings on coronaviruses, together 
with data on henipavimses (23-25, 28), sug¬ 
gest that genetic diversity exists among zoo¬ 
notic viruses in bats, increasing the possibility 
of variants crossing the species barrier and 
causing outbreaks of disease in human popula¬ 
tions. It is therefore essential that we en¬ 
hance our knowledge and understanding of 
reservoir host distribution, animal-animal and 
human-animal interaction (particularly within 
the wet-market system), and the genetic diversity 
of bat-bome viruses to prevent future outbreaks. 
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Hypothalamus of Adult Mice: 
Potential Role in Energy Balance 
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Ciliary neurotrophic factor (CNTF) induces weight loss in obese rodents and 
humans, and for reasons that are not understood, its effects persist after the 
cessation of treatment. Here we demonstrate that centrally administered CNTF 
induces cell proliferation in feeding centers of the murine hypothalamus. Many 
of the newborn cells express neuronal markers and show functional phenotypes 
relevant for energy-balance control, including a capacity for leptin-induced phos¬ 
phorylation of signal transducer and activator of transcription 3 (STAT3). Co¬ 
administration of the mitotic blocker cytosine-fi-D-arabinofuranoside (Ara-C) 
eliminates the proliferation of neural cells and abrogates the long-term, but not 
the short-term, effect of CNTF on body weight. These findings link the sustained 
effect of CNTF on energy balance to hypothalamic neurogenesis and suggest that 
regulated hypothalamic neurogenesis in adult mice may play a previously 
unappreciated role in physiology and disease. 


The obesity epidemic has prompted major ef¬ 
forts to develop safe and effective therapies 
(1, 2). However, approved drugs for obesity 
have limited efficacy and act only acutely, 
with patients rapidly regaining weight after 
terminating treatment (3). Only the neuro¬ 
cytokine ciliary neurotrophic factor (CNTF) 
and Axokine, an analog of CNTF developed 
as a drug candidate for the treatment of obe¬ 
sity, appear to deviate from this paradigm. 
Rodents and patients treated with Axokine 
were reported to maintain lowered body 
weights weeks to months after the cessation 
of treatment (4, 5). This feature of Axokine/ 
CNTF action is unexplained and suggests that 
CNTF induces long-lasting changes in one or 
more elements of the energy-balance circuitry. 

In rodents, CNTF is most potent when 
administered directly into the cerebrospinal 
fluid (d) and activates signaling cascades in 
hypothalamic nuclei involved in feeding 
control (5, 7, 8). For instance, CNTF acti¬ 
vates phosphorylation of signal transducer 
and activator of transcription 3 (STAT3) in a 
population of hypothalamic neurons that 
substantially overlaps with those activated 
by leptin (5). However, in contrast to CNTF, 
leptin-treated animals do not maintain their 
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lowered body weight after the cessation of 
treatment. We thus sought a CNTF-specific 
mechanism to explain this long-term effect. 

CNTF supports the survival of neurons in 
vitro and in vivo (9) and has also been impli¬ 
cated in the maintenance of adult neural stem 
cells (10). Furthennore, other trophic factors, 
such as epidermal growth factor and fibroblast 
growth factor 2, are known to act as mitogens 
on adult neuronal progenitors (11, 12), and they 
promote the functional regeneration of hippo¬ 
campal pyramidal neurons (13). Neurogenesis 
in the adult brain is most clearly defined in the 
subventricular zone (SVZ) of the lateral ven¬ 
tricles and the subgranular zone (SGZ) of the 
hippocampal fonnation (14). However, recent 
reports indicate that the neuroproliferative 
potency in the adult extends to other brain 
structures, including the hypothalamus (15-17). 
On the basis of these findings, we hypothesized 
that the long-term effect of CNTF on body- 
weight regulation might involve neurogenesis 
in the hypothalamus, which is the brain region 
most relevant for energy-balance regulation. 

To assess the mitogenic potency of CNTF 
in the adult nervous system in vivo, we 
delivered the cell-proliferation marker bromo- 
deoxyuridine (BrdU) alone (vehicle treatment) 
or together with CNTF directly into the 
cerebrospinal fluid of mouse brains (18). 
CNTF and BrdU were continuously infused 
for 7 days into the right lateral ventricle using 
osmotic minipumps. Mice were switched to a 
high-fat diet two months before surgery and 
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were kept on this diet throughout the experi¬ 
ments. In accordance with previous results (5), 
CNTF-treated mice showed a marked reduc¬ 
tion in body weights (Fig. 1A), which persisted 
after termination of CNTF delivery. Mice were 
killed 22 days after surgery, and brain sections 
were immunostained with an antibody against 
BrdU. Because BrdU incorporates into DNA of 
dividing cells, BrdU-positive (BrdU + ) cells are 
thought to represent newborn cells. Figure IB 
shows coronal sections of vehicle- and CNTF- 
infused animals at the level of the arcuate, 
ventromedial, and dorsomedial nuclei, well- 
known hypothalamic centers for energy-balance 
regulation (19). In vehicle-infused animals, few 
BrdU + cells were detected in the parenchyma 
surrounding the third ventricle (Fig. IB, left). 
Administration with CNTF led to a dramatic 
increase of BrdU+ cells (Fig. IB, right). Note 
the higher density of BrdU + cells at the base 
of the third ventricle, which is part of the ar¬ 
cuate nucleus/median eminence. 

The pattern of CNTF receptor (CNTFR) 
mRNA expression is consistent with this ob¬ 
servation. In situ hybridization using a ribo- 
probe against CNTFR mRNA revealed strong 
staining in the walls of the basal third ventricle 
and surrounding arcuate nucleus parenchyma 
(Fig. 1C). Because this section originated from 
an animal treated with both CNTF and BrdU, 
we colabeled with antibodies to BrdU. Many 
BrdU + cells were positive for CNTFR ex¬ 
pression, indicating that CNTF, at least in part, 
directly promotes cell division by binding to 
CNTFR on putative neural progenitor cells 
(Fig. ID, inset). By counting all newly gener¬ 
ated cells in the caudal hypothalamus, CNTF 
treatment led to a marked increase of BrdU + 
cells over vehicle-infused animals (Fig. IE). 
The total number of BrdU+ cells in CNTF- 
treated animals remained constant for at least 
2 weeks after the infusion period. Subsequent¬ 
ly, the numbers decreased but plateaued at a 
high level. Vehicle-infused animals showed a 
similar fractional decrease over time. Thus it 
appears that the majority of hypothalamic 
BrdU + cells do not die or migrate to distant 
areas as reported for newborn neurons of the 
SVZ, which follow the rostral migratory stream 
toward the olfactory bulb (20). 

To investigate the origin of adult-bom cells 
in the hypothalamus, we examined CNTF and 
vehicle-infused brains every 12 hours starting 
48 hours after surgery, a time when the in¬ 
fused CNTF/BrdU should just reach the ven¬ 
tricular system (18). Hypothalamic BrdU 
incorporation was first detected 60 to 72 hours 
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